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Weigh the benefits 
of fuzzy-logic vs classical 
control in a disk-drive spindle 



Brian P Tremaine, Seagate Technology 



You can apply fuzzy- and classical-control 
techniques to any servo-control loop. Which 
technique you use depends heavily on the 
nonlinearities in a system. 

| A furor is raging over whether you should use fuzzy-logic or 
classical-control tools to analyze a complex servo-control sys- 
tem Fuzzy-logic advocates claim it eliminates the need for a 
mathematical system description, while classical-control 
advocates claim fuzzy logic's lack of analysis tools makes it 
undesirable. To observe these extreme views, compare the 
classical- and fuzzy-control analysis of a disk-drive dc spin- 
dle motor. The comparison uses a C program, which simu- 
lates the plant and classical feedback-control signal 
The FIDE (fuzzy inference-development) program from 

Aptromx Inc compiles the plant program to simulate a fuzzy 

controller. The spindle-motor example 

shows that even a basic system can have f ~ ~ 

nonlinearities, which preclude linear 

analysis. However, knowledge of the 

system equations can aid the design of 

a fuzzy controller. 
The dc spindle motor, driver elec- 

■ tronics. and digital transducers com- 

■ prise the system's plant and are inher- 

■ tntly nonlinear. The driver develops 

■ only positive torque, which is a common 
■practice in disk-drive applications. 
■Therefore, the motor depends only on 
■coulomb friction to decelerate. The 
■plant's nonlinearities include satura- 
tion of the error signal, torque control 
■in only one quadrant, and quantization 
■of the transducer signals from the out- 
■put control voltage. Therefore, linear- 
■analysis methods do not apply for either 
■the classical or fuzzy set of control 
■tools. This analysis compares the dif- 
ferences between a classical propor- 
tional plus integral (PI) controller and a 
■uzzy controller. 



The plant model (Fig 1) is a dc 3-phase, 12-pole spindle 
motor along with the driver electronics. The output circuit- 

hTT™ three half - brid e e drivers with current-sense 
feedback. The input to the driver is an analog voltage which 
generates the spindle motor current through a current 
amplifier haying a gain of gm (amps/volt). When the input 
voltage equals some reference voltage, the current command 
is zero and the output stages are off. During this time 
coulomb friction decelerates the motor. Maximum current 
drive occurs when the input voltage is at V. The motor resis- 
tance, supply voltage, and the back EMF from the motor 
limit the maximum current drive to the motor, 

Simniator doesn't model current control 

The current-control loop around the motor and the current 
amplifier typically have a bandwidth several orders of mag- 
nitude greater than the motors velocity bandwidth Thert 
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fore, che plant simulation doesn't include the dynamics of the 
current-control loop in Fig 1. Instead, the plant model calcu- 
lates the maximum available current and drives the motor 
using a current source, which is equal to either the com- 
manded current or the maximum available currerft, whichev- 
er is smaller. 

The classical-control analysis applies linearization to the 
nonlinear plant to design a proportional PI controller. The 
linear model does not predict the correct simulated tran- 
sient response. The simulated plant model calculates 
motor speed in rad/sec, position in rad, and a digital 
encoder pulse, called f . A digital tachometer and digital 
phase detector use the f com signal, which occurs nine 
times/revolution of the motor (N p =9). The linearized 
transfer function of the motor and current driver in con- 
tinuous time is: 



H(s)=(gm-Kt/Jm)/(s + Kv/Jm), rad/sec-V, 



CD 



where gm=0.5 A/V. transconductance; Kt=1.05 oz-in./A, 
torque constant; Kv=0.001 oz-in./sec, viscous constant; and 
Jm=752E-6 oz-in.-sec-, inertia. 

The equation represents the transfer function from the 
input voltage to the motor's radial velocity output. 
Because the tachometer's outputs are discrete, the con- 
tinuous-time representation is converted to discrete time 
using the z-transt'orm. The transform involves an approx- 
imation because the sample interval is at a fixed motor 
rotation (2-n7N p ) and not at a fixed sample time. As a prac- 
tical matter, the discrete system is often analyzed assum- 
ing T is a fixed sample time. Although we make this 
approximation, linearizing the system and assuming a 
fixed sample time completely ignores the real system's 
nonlinearity. Stability of the linearized system does not 
ensure global stability of the real system. An appropriate 
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Fig 2 — Driving a spindle motor from o dead stop to a final spindle speed using o eli 
proportional integral (PI) controller shows a severe target overshoot. The integrator 
off the motor current at Paint 2 to allow the motor to settle using coulomb friction. 
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Lyaponov function could determine a nonlinear systei 
degree of stability. Finding such a Lyaponov function 
difficult, however. 

You can derive a difference equation directly from Eq 1 
assuming a small time increment, as follows: 

Mt+dt)=t l-dt(Kv/Jm)).<o(t)+(gm.Ktdt/Jm)-v(t). (2 

If dt=T and t=kT. then Eq 2 becomes 

(iu(kT+T=(l-T(Kv/Jm))-(u(kT)+(gm-KtT/Jm).v(kT). 

Noting that the z-transform of uj(kT+T is u>(z)-z, the oven 
transfer function is 

H(z)=gm-Kt T/Jm)/(z-(.l-T-Kv/Jm)), rad/sec-V. 

Digital tachometer provides feedback 

Regardless of whether the controller is classical or 
you must measure the spindle speed for feedback control. Fij 
1 includes a block diagram labeled "digital tachometer"showfl 
ing the measurement of the spindle speed, uCt), The outputs 
of this block is xd_err, a digital word representing speed! 
error. 

The digital tachometer counts the period of the f^ pulse* 
with 1-iisec resolution (t ( ). You subtract the final count from! 
a reference count of 1880. The tachometer then truncates thel 
count to fall within -128 to +127 counts to generate thel 
xd_err signal. If the time interval is T t and the average veloc* 
ity over an t m period is w then wJTJ is the average? 
velocity over the time interval: 

»„(T,)-2^/(N„T t ). 

This equation is nonlinear in the variable T.. To linearize I 
equation, you can write a Taylor series around the referenc 
velocity to^and the reference period T t . You linearize a func-i 
tion using a Taylor series about an equi- 
librium point. If you linearize f(x,t) 
about the point xO. then the first 
terms of the Taylor series are 

iXx,t)=ftxO,t)+f'(x0,tXx-x0), 

where f denotes the derivative witl 
respect to x. Applying this to Eq S 
yields 

<«„ < lT> u „r(„ i yrj'<T,-T„,). (7)1 

Because w^T^ is the average velocity; 
over the period T„ you can express this 
equation in terms of the velocity at thi 
beginning and the end of the period, a 
follows: 

w„,(k)=(l/2)-( w <k) + ( 1 >(k-l)). 

Combining Eqs 7 and 8 yields 

(colki-uj^+twtk-D-w^)- 

Denoting the difference between w and ; 
tu^ as ui, and measuring the time interH 
val in terms of t reduces 
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We've always been known for our high-quality wiremount connectors. Shouldn't the 
world know about our boardmount connectors, too? 

Well, we can get you a complete line of boardmount products: shrouded headers, 
in both low-profile and latch-eject models; shunts; pin strips. Polarized sockets, straight 
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We can get you high- temperature bodies, true surface-mounts and through-boards, 
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<ue(k)+ 0) e<k-l=-2 (i« ri /rj t | .d(k), (10) 

where d(k) is the digital count of the time interval. 

The z-translbrra for Eq 10 is the transfer function of the 
digital tachometer, as follows: 

me(z)+ < ae(z)/z=-2((« i yT n( )(T r t l -d(z). 
D(z)=d(z)/o«(z)=-T m /(2-»^. t,)). 1 +z)/z. 

If the motor velocity is signilicantly different from the refer- 
ence frequency (u^), the Taylor series expansion is no longer 
accurate. In addition, if the velocity is low, the sample time 
is long, and the phase delay of D(z) is long. 

The plant-transfer function in Eq 1 has a maximum phase 
delay of 90". Even when you represent Eq 1 in discrete form 
(Eq 4), the closed-loop bandwidth is sufficiently below the 
Nyquist frequency, 1/T, that the phase delay caused by sam- 
pling is small. Therefore, you can stabilize the plant using 
only proportional feedback from the tachometer. However, 
the speed error is not exactly zero. Offsets in the driver stage 
and motor running torque relative to the driver input 
requires a speed error to counteract the offsets. 

To force the speed error to zero, classical feedback-control 
systems sum the speed error and the time integral of the 
speed error (PI) to generate the feedback-control law. The 
control-law calculation for a discrete time system follows: 

Vint(k)=Vint(k-l*k,.xd_err(k); (13) 



The open-loop gain of the linearized system is 

H(z).C(z)-D(z), Gol(z)=D(z).((k„+k.).gm K. T/J > 
(z-k [ )/(k t 4.k,))/[z.(l-T.Kv/Jm)).(z-l)]. ' " (i 7) 

Substituting some of the system parameters into Eq 17 
yields open-loop poles as zl=l, and z2=0.99. The digital 
tachometer produces a pole at 8-0 and a zero at z=-l. The 
ratio of the integral gain to proportional gain determines the 
zero at zsl/d+k/k^. 

The linearized model predicts closed-loop poles that are 
well-damped for small k/k p and under-damped for large k/k . 
This result is in contrast to the simulation results (Fig 2) and 
highlights the fact that nonlinearities in a system can inval- 
idate the use of linear methods. 

Fig 2 shows a plot of the simulated spindle speed using a 



Voutftlsk^xd.errM+VinMk), 



(14) 

where xd_err is the digital word representing speed error. 
Eq 13 is the difference equation for the digital integrator, 
and Eq 14 is the linear sum of the speed error and time inte- 
gral of the speed error. Y ou can combine Eqs 13 and 14 into 
a transfer function in the z-domain for the equivalent digital 
controller, C(z), as follows: 

C(z)=Vout(z)/xd_err(z)=kp+ki-z/(z-l). 
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3 SPINDLE disk drive ipindJc PI control 

J bn» on™™, 6 July 1993. MM 1 Isa 1994 

>vi_or-bio-:-l2»0127 [ 

Ntfjjjrg. (8-128.255. 9-64.0), 
*e_Med (9-128.0. 9-64,235. 90.0). 
Zoo (9.64.0. 90.255. 964,0), 

Pos_l*d (90.0. 964J55. 9127.0). 
Po._L»ie (964.0. 9127.253)]; 

>v_r*l-bili-fl0253 ( 

NesJUrjj. (90,255. 964,0). 

Ne^Med (90.0. 964.255. 9128.0). 
Zoo (964,0, 9 128 .255. 9192,0). 

Pr*_rVted (9128,0. 4192.255 . 9233,0). 

Pm_L > rgo(9l92.0. 9255.255)]; 

<aror-biB-: 00235«( 

0J>. 
I 

Zero- 128, 
Pr*_Smiil - 248. 

■ 255); 



<v_new "bio:-: 00255' ( 
Nef_L"rse. 0.0. 
Noj_M«l . 64. 
Zoo- 128. 
Pos.Med. 192, 
Prajjoje . 255 ); 

t 

if xd_err is Ne g_L*ree 
ifwJ_eniiNeg_M«i 

if rd_crr is Zero and v_old is Poa Large | 
if xd_crr is Zero arsl v jjtd la Pos^Med I 
if rd_err is Zero and v old is 7cn> 
if rdjsr is Zero and v old is Neg_Med t 
if id erris Zero and v_old is Ncg_Laree t 
ifad.err if Pw_Mod raid v.old is Poa Med t 
if id_or is Poa_Med and v.old ii Zoo 
if xd_or is Pos Jflod and v.old is StfJM 
if sd.etr is Pos.Med and v.old is Nojjjige 
ifxd.errisPosJjrgc 
3 

if sd_err is Nef_Lsrge 
ifsd_errisNe8_Med 
ifxd.emsZero 
if wi_ejT is Poa_Med 
if sd^rar is Pos J-argc 



ov_newisPosJjrge; 
ov_iirr»«r^._Med; 
M vjiew is Zoo: 
n v_new is Nef^Med; 
a v_aew is rVeR_Lsrge; 
n v_ocw la Pos_Large; 
» v_new is Pos_MetL 
■ vjsaw is Zero: 
sa vjsew is (seized; 
!6rav_i>t»iaPcsJjujes 

□>=«oro>UNe t _Ur8o; 

daeo error is Net>fad; 

daeo error is Zero; 
(fien error is Pos_Small; 
Ifaao error is Pos_Large; 
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Z'7Lt^lZi° Proportioniil pi.. m.egrol »»- ,3 ; ampiled with Aplroni,', FIDE (hiiiy inference development) 
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fig 5—~The rule table for the fuzzy integral controller contains redundant rows of rules that 
can be condensed into fewer rules for the fuzzy-inference unit to follow. 
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analogous to the proportional plus inte- 
gral control block of the classical 
method. Because an FIU has no memo- 
ry, fuzzy control accomplishes the inte- 
gral function by using the feedback 
from the FIU's output to its input. The 
FIU has the FIDE listing shown in Fig 
4. The inputs to the FIU are the speed 
error (xd_err) and previous output 
integral count (V_old). The outputs of 
the FIU are proportional count (error) 
and integral count (V_new). Fig 5 
shows a table representation of the FIU 
rules. The control law is a summation of 
both outputs 

Vout=16-(error-128)/255+3 
(v,new-l28)/255+3. 

Fig 6 shows a plot of motor speed and 
current for spin-up from a dead stop. 
The steady state speed nearly equals 
the commanded speed, as in the inte- 
gral controller. However, the fuzzy 
method improves transient response, 
decreases overshoot, and provides 
faster settling. A classical PI system 
involves a tradeoff between overshoot 
and settling time. The fuzzy-system 
rules and membership functions effec- 
tively implement "anti-wmdup" inte- 
gral control. When xd_err is positive, 
(over speed) the integral term (v_new) 
is forced positive to zero the spindle- 
motor current immediately to cause 
deceleration. 

Table 1 details the fuzzy model: 



fig 6— The spin-up of the spindle motor from a dead stop using fuzzy integral control shows 

|.„ au>»L>^* iL._ ,..;iL m Dl ■ ■■■■■ U. ■ tLji «Ju»> 

less oversnooi man wtin a classical ri controller, me integrator reduces motor current io 
zero when the system is over speed to cause deceleration. 
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classical PI controller with a nonlinear plant. In Eqs 13 and 
14, kp.3/2048, and k p =48/2048. The speed overshoots the tar- 
get and slowly decays with a linear ramp to the target speed 
value. From time t=0 until Point 1. the speed error is nega- 
tive, and the integrator continues to accumulate until it sat- 
urates. Not until the speed error begins to go positive does 
the integrator begin to come out of saturation. At Point 2 the 
integrator has recovered enough to allow the control voltage 
to command zero current. From Point 2 until Point 3 the only 
deceleration is due to coulomb friction. Beyond Point 3, the 
system settles to a steady state value. 

Fuzzy control uses different control blocks 

Fuzzy integral control analyzes the spindle motor using a 
fuzzy-inference unit (FIU) (Fig 3). This analysis block is 
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The execution time for the above model with a 68HC11 
M.P is about 1.25 msec using a 50-nsec bus cycle time. The 
execution time is sufficient for the sample time and f cnB 
period of 1.88 msec. The model is optimized for the mini- 
mum number of rules, which is an important step affect- 
ing the choice of hardware. The initial pass of the FIU 
has 30 rules and 877 code bytes. Although this pass 
executes in 1.6 msec, there is little margin in processor 
bandwidth. 

The output v_new is a function of xd.err and v_old. 
Because each input has five labels, you can define v_new as 
a 5x5 matrix, or 25 rules. The output, error, is a function 
of xd_err only. Because xd_err has five labels, five rules 
can define the output, error. Therefore, there are 30 rules 
for the the first pass. In the last row of Fig 5, when xd_err 




UPS Standard 



UPS Expedited 



UPS Express 



How small a continent do you want to make it? 



^uingiy. ror more information on how we can rsiW M„ ,k a ■ , 1 1 
™^«"»a.l-800-PICKH^ 

CIRCLE NO. 163 

H July 7, 1994 ■ 143 



EDN-Design Feature 



FUZZY AND CLASSICAL CONTROL 

is P08_Large, V_new is Poa_Large, regardless of v_old. 
This condition lets you replace five rules with one rule: "If 
xd.err is Pos.Large, then v_new is Pos_Large." Modifica- 
tions such as these can reduce the number of rules from 30 
to 17. , 

In comparing a classical PI servo-control loop with a 
fuzzy-control loop for a dc spindle motor, the plant and the 
tachometer have inherent nonlinearities that preclude 
applying linear stability theory, regardless of the control 
method. Knowledge of classical PI servo-control principles 
helps with the design of a fuzzy controller. The perfor- 
mance in terms of the spin-up settling time and the tran- 
sient response is superior using the fuzzy system compared 
with the classical PI controller when there are severe non- 
linearities. 

The C source code for both the FIDE fuzzy modules and 
the plant model is available from the author. In addition to 
the PI controller, you can also use the Aptronix bulletin- 
board system to access source and executable simulation 
code for proportional control, PLL control, and dual-phase 
and tachometer loops for classical and fuzzy control. Wkw 



2. Franklin, G, D Powel, and M Workman, Digital Con- 
trol of Dynamic Systems, Second Edition, Addison Wesley, 
1990. 
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